Decoherence during Inflation: the generation of classical inhomogeneities 
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We show how the quantum to classical transition of the cosmological fluctuations produced during 
inflation can be described by means of the influence functional and the master equation. We split 
the inflaton field into the system-field (long- wavelength modes), and the environment, represented 
by its own short-wavelength modes. We compute the decoherence times for the system-field modes 
and compare them with the other time scales of the model. We present the renormalized stochas- 
tic Langevin equation for an homogeneous system-field and then we analyze the influence of the 
environment on the power spectrum for some modes in the system. 
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I. INTRODUCTION 

The emergence of classical physics from quantum be- 
haviour is important for several physical phenomena in 
the early Universe. This is beyond the fundamental re- 
quirement that only after the Planck time can the metric 
of the Universe be assumed to be classical. For example, 
the inflationary era is assumed to have been induced by 
scalar inflaton fields, with simple potentials Such 
fields are typically assumed to have classical behaviour, 
although in principle a full quantum description should 
be used. In fact, the origin of large scale structure in 
the Universe can be traced back to quantum fluctuations 
that, after crossing the Hubble radius, were frozen and 
became classical, stochastic, inhomogeneities ||. 

It is generally assumed that several phase transitions 
have occurred during the expansion of the Universe £| . 
As in the case for the inflaton fields, the (scalar) or- 
der parameter fields that describe these transitions are 
described classically. However, the description of early 
universe phase transitions from first principles is intrin- 
sically quantum mechanical . As a specific application 
[f| of the previous point, the very notion of topological 
defects (e.g. strings and monopolcs) that characterize 
the domain structure after a finite-time transition, and 
whose presence has consequences for the early universe, 
is based on this assumption of classical behaviour for the 
order parameter Q, as it distributes itself between the 
several degenerate ground states of the ordered system. 

In previous publications, one of us has analyzed the 
emergence of a classical order parameter during a second 
order phase transition and the role of decoherence in the 
process of topological defect formation 0, |jj 0, |n| . 

In the present paper our concern is directly related 
with the first point above, the quantum to classical tran- 
sition of the inflaton. Any approach must take into ac- 
count both the quantum nature of the scalar field and the 



non-equilibrium aspects of the process [12 ] . The problem 
of the quantum to classical transition in the context of 
inflationary models was first addressed by Guth and Pi 
|13|. In that work, the authors used an inverted harmonic 
oscillator as a toy model to describe the early time evo- 
lution of the inflaton, starting from a Gaussian quantum 
state centered on the maximum of the potential. They 
subsequently showed that, according to Schrodinger's 
equation, the initial wave packet maintains its Gaussian 
shape (due to the linearity of the model) . Since the wave 
function is Gaussian, the Wigner function is positive for 
all times. Moreover, it peaks on the classical trajectories 
in phase space as the wave function spreads. The Wigner 
function can then be interpreted as a classical probability 
distribution for coordinates and momenta, showing sharp 
classical correlations at long times. In other words, the 
initial Gaussian state becomes highly squeezed and indis- 
tinguishable from a classical stochastic process. In this 
sense, one recovers a classical evolution of the inflaton 
rolling down the hill. 

A similar approach has been used by many authors 
to describe the appearance of classical inhomogeneities 
from quantum fluctuations in the inflationary era [I3 . ll5| . 
Indeed, a massless free field <fi in an expanding universe 
can be written as <fi — a~ lr ip where a is the scale factor 
and the Fourier modes of the field ip satisfy the linear 
equation 



(k 2 = o. 

a 



(1) 
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For sufficiently long- wavelengths (k 2 <C a" /a), this equa- 
tion describes an unstable oscillator. If one considers an 
initial Gaussian wave function, it will remain Gaussian 
for all times, and it will spread with time. As with the 
toy model of Guth and Pi, one can show that classical 
correlations do appear, and that the Wigner function can 
again be interpreted as a classical probability distribution 
in phase space. (It is interesting to note that a similar 
mechanism can be invoked to explain the origin of a clas- 
sical, cosmological magnetic field from amplification of 
quantum fluctuations). 

However, classical correlations are only one aspect of 
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classical behaviour. It was subsequently recognized that, 
in order to have a complete classical limit, the role of the 
environment is crucial, since its interaction with the sys- 
tem distinguishes the field basis as the pointer basis [l6| • 
[We are reminded that, even for the fundamental prob- 
lem of the space-time metric becoming classical, simple 
arguments based on minisuperspace models suggest that 
the classical treatment is only correct because of the in- 
teraction of the metric with other quantum degrees of 
freedom [HI-] 

While these linear instabilities cited above characterize 
free fields, the approach fails when interactions are taken 
into account. Indeed, as shown again in simple quantum 
mechanical models (e.g. the anharmonic inverted oscil- 
lator), an initially Gaussian wave function becomes non- 
Gaussian when evolved numerically with the Schrodinger 
equation. The Wigner function now develops negative 
parts, and its interpretation as a classical probability 
breaks down [lOj. One can always force the Gaussianity 
of the wave function by using a Gaussian variational wave 
function as an approximate solution of the Schrodinger 
equation, but this approximation deviates significantly 
from the exact solution as the wave function probes the 
non-linearities of the potential [ToL Il8| . 

When interactions are taken into account, classical be- 
haviour is recovered only for "open systems", in which 
the observable degrees of freedom interact with their en- 
vironment. When this interaction produces both a diag- 
onalization of the reduced density matrix and a positive 
Wigner function, the quantum to classical transition is 
completed [T^ |. 

In Ref. |10| has been studied an anharmonic inverted 
oscillator coupled to a high temperature environment. 
Under some considerations, it was shown that the system 
becomes classical very quickly, even before the wave func- 
tion probes the non-linearities of the potential. Being an 
early time event, the quantum to classical transition can 
now be studied perturbatively. In general, recoherence 
effects are not expected 20|. Taking these facts into ac- 
count, we have extended the approach to field theory 
models H, 0- ^ n ne ld theory, one is usually interested 
in the long-wavelengths of the order parameter. Even 
the early universe is replete with fields of all sorts which 
comprise a rich environment, in the inflationary example, 
we considered a model in which the system-field interacts 
with the environment-field, including only its own short- 
wavelengths. This is enough during inflation. Assuming 
weak self-coupling constant (a flat inflaton potential) we 
have shown that decoherence is an event shorter than the 
time i C nd, which is a typical time-scale for the duration 
of inflation. As a result, perturbative calculations are 
justified 0. Subsequent dynamics can be described by 
a stochastic Langevin equation, the details of which are 
only known for early times |2l| . 

In our approach, the quantum to classical transition 
is defined by the diagonalization of the reduced density 
matrix. In phase transitions the separation between long 
and short-wavelengths is determined by their stability, 



which depends on the parameters of the potential. Dur- 
ing Inflation, this separation is set by the existence of the 
Hubble radius. Modes cross the apparent horizon during 
their evolution, and they are usually treated as classical. 
The main motivation of this paper is to present a for- 
mal way to understand this statement within the open 
quantum system approach. In the last sense, decoher- 
ence is the critical ingredient if we are to dynamically 
demonstrate the quantum-to-classical transition of the 
open system. 

The splitting between short and long-wavelength 
modes may be done using a time-dependent or time- 
independent comoving cut-off as well as a smoother win- 
dow function [2lll22j . Since any time-dependent splitting 
produce an effective and arbitrary (split-dependent) in- 
teraction between the system and environment degrees 
of freedom, which can only be discarded by making some 
additional assumptions, we consider convenient to use a 
time-independent comoving cut-off, as it has been done 
in Refs. [$1 |3 

The paper is organized as follows. In Section II we 
introduce our model. This is a theory containing a real 
system-field, massless and minimally coupled to a fixed 
de Sitter background. We compute the influence func- 
tional by integrating out the environmental sector of the 
field, composed by the short-wavelength modes. Section 
III is dedicated to reviewing the evaluation of the master 
equation and the diffusion coefficients which are relevant 
in order to study decoherence. In Section IV we analyze 
the diffusion coefficients and evaluate upper bounds on 
the decoherence times. As we will see, decoherence takes 
place before the end of the inflationary period. Section V 
is concerned with the effective stochastic evolution of the 
system. We present the renormalized stochastic Langevin 
equation for an homogeneous system-field and then we 
analyze the influence of the environment on the power 
spectrum for some modes in the system. Section VI con- 
tains our final remarks. Two short appendices fill in some 
of the detail. Throughout the paper we use units such 
that h = c = 1 . 



II. THE INFLUENCE FUNCTIONAL AND THE 
DENSITY MATRIX 



Let us consider a massless quantum scalar field, mini- 
mally coupled to a de Sitter spacetime ds 2 = a(r])[dri 2 — 
dx 2 } (where r\ is the conformal time, di] = dt/a(t) with 
t the cosmic time), with a quartic self-interaction. The 
classical action is given by 



sr. 



d 4 x a 4 (?y) 



2a 2 (rj) 2a 2 (77) 



\<j) 



(2) 



where a(ry) = —l/(Hr]) and = d(j)/dr] {a{r]i) — 1 [77,; = 
—H -1 ], and H' 1 is the Hubble radius). Let us make a 
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system-environment field splitting [24 



The influence functional (or Feynman- Vernon functional) 



is defined as 



(3) 



where the system-field (f> < contains the modes with wave 
vectors shorter than a critical value A = 2tt/X c , while the 
environment-field 0> contains wave vectors longer than 
A. As we set a{rji) — I , a physical length A p h ys = a(rj)X 
coincides with the corresponding comoving length A at 
the initial time. Therefore, the splitting between system 
and environment gives a system sector constituted by all 
the modes with physical wavelengths shorter than the 
critical length A c at the initial time r\i. 

After splitting, the total action can be written as 



F\ 



P>\ 



>>i'<P>i>Vi. 



S[</>] - s< 



s >[ 



5';, 



(4) 



where Sq denotes the free field action and the interaction 
terms are given by 

Stat[<t><,<h] = -A / d 4 x a 4 ( V ) {<j>\{x) + 4 (x) 
+ 60 2 < {x)<t>l (x) + 403 (x)</>> (x) + (x)4>l (x) } .(5) 

The total density matrix elements (for the system and 
environment fields) are defined as 



Vcj>+ / Vr> exp {i(S[c/>+] + S int [</>+, </>+])} 
J( t>>i 

x exp{-*(S[0>] + Stat^-, 0;])} . (fl) 

This functional takes into account the effect of the en- 
vironment on the system. The influence functional de- 
scribes the averaged effect of the environmental degrees 
of freedom on the system degrees of freedom to which 
they are coupled. With this functional, one can identify 
a noise and dissipation kernel related by some kind of 
fluctuation-dissipation relation. This relation is impor- 
tant when one is interested in possible stationary states 
where a balance is eventually reached. During inflation 
we have the field (inflaton) on a very flat potential, away 
from its minimum, and we are, in general, only interested 
in the dynamics over some relatively small time. For ex- 
ample, we would neglect dissipation during the slow-roll 
period; but it is not correct during the eventual reheating 
phase. 

We define the influence action S A \4><, </><] and the 



|0<, 4> > ;v\ = (<M>IpMI<M>), ( 6 ) coarse grained effective action (CGEA) A[</>^,^] as 



where \<j><) and \<f>>) are the eigenstates of the field op- 
erators <j) < and </>>, respectively. For simplicity, we will 
assume that the interaction is turned on at the initial 
time rji and that, at this time, the system and the envi- 
ronment are not correlated (we ignore, for the moment, 
the physical consequences of such a choice, it has been 
discussed in Therefore, the total density operator 

can be written as the product of the density operator for 
the system and for the environment 



F\ 



;] = exp{iSA[4>^, 



A[tf>+, <t><] = S[<t>+] - S[r<] + SA[<p+, 



(12) 



(13) 



p[m] = p>[m]p<[m\- 



(7) 



We will further assume that the initial state of the envi- 
ronment is the Bunch-Davies vacuum pl|. 

We are interested in the influence of the environment 
on the evolution of the system. Therefore the reduced 
density matrix is the object of relevance. It is defined by 

pA<f>t\<f><;v\= /^>p[^,0>|0<,^>;»7]. (8) 

The reduced density matrix evolves in time by means of 
##</#</;»?] = J J p r [4>%i\4><ilVi] (9) 

where J r is the reduced evolution operator 

r4>t f r*< f 



We will calculate the influence action perturbatively in A 
and we will consider only terms up to order A 2 and one 
loop in the h expansion. The influence action has the 
following form |24|: 

0<] = (5 int [0+,0+]) o -(5 int [^,0<]) o (14) 

- «(S'int[0>,0<]<5mt['/'>,0<])o 
+ i(S'int[0>,0<])o(5'int[0>,0<])o 

+ ^{(5 int [0+,0+]) 2 -(5 int [0+,0+] 2 ) o } 

+ 2 {(' S wt[0>J0<])o - (Sint [4>>, 0<] 2 )o} 

where ( )o is the quantum average with respect to the 
free field action of the environment, defined as 



(15) 



Dj>+ / D4>-B[^^-} 

x exp {i (s> [</>+} - s>[<t>-})} (4>Up>hU; 



x e X p{i(S<[0+]-S<[0<])}F[0+ </>-]. (10) 



Here p> is the Bunch-Davies vacuum state assumed for 
the environment. 
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The influence functional can be computed, and the re- 
sult is 

Re5A = -X J d 4 Xl a 4 (7y) {2A 4 (iei) 

- 12A 2 (x 1 )iGl + (x 1 ,x 1 )} 

+ A 2 J d A xi J d i x 2 a 4 (r]i) a i (i] 2 ) 6(771 ~ m) 

x {64A 3 (xi)ReG^ + {xi,x 2 )Z 3 (x 2 ) 

+ 288A 2 (x 1 )ImGf + {x 1 ,x 2 )^2{x 2 )}, (16) 

haSA = X 2 J d 4 Xl J d 4 x 2 a 4 (?7i) a 4 (?7 2 ) 

x {32A 3 (x 1 )ImG^ + (a;i,x 2 )A3(x 2 ) 

144A 2 (a: 1 )ReGf + ( a;i ,a ;2 )A 2 (a ;2 )}, (17) 

where Xj denotes (r)j,Xj), Q(x) is the Heaviside step 
function, and the integrations in time run from rji to rj. 
G+ + (xi,x 2 ) = i(0>(xi)0>(a: 2 ))o is the relevant short- 
wavelength closed time-path correlator (it is proportional 
to the Feynmann propagator of the environment field, 
where the integration over momenta is restricted by the 
presence of the infrared cut-off A), and we have defined 



+ (18) 



with n = 1, 2, 3. 



III. MASTER EQUATION AND DIFFUSION 
COEFFICIENTS 

In this Section we obtain the evolution equation for 
the reduced density matrix (master equation), paying 
particular attention to the diffusion terms, which are re- 
sponsible for decoherence. To do so, we closel y fo llow 
the quantum Brownian motion (QBM) example |2al27j| . 
translated into quantum field theory 0, 0] ■ 

The first step in the evaluation of the master equa- 
tion is the calculation of the density matrix propagator 
J r from Ea. l|10[l . In order to solve the functional integra- 
tion which defines the reduced propagator, we perform a 
saddle point approximation, assuming the classical field 
configuration dominates functional integrals, 



(19) 



where <^< c i is the solution of the semiclassical equation of 



motion 5ReA/6<j)<\ 



= with boundary conditions 



0< c l(»7t) = (f>%i an d ^div) = </></• Since we are working 
up to A 2 order, we can evaluate the influence functional 
using the solutions of the free field equations. This classi- 
cal equation is <f) < +2W</> < — V 2 </>< = 0, (H = a' '(rj) / a(rj)) . 
A Fourier mode ipz of the field ijj = a{rj)4> <1 satisfies 



k 2 --)^, = 0, 



(20) 



where we have used the fact that a" /a = 2/rf. It 
is important to note that for long-wavelength modes, 
k <C 2/ry 2 , Eq. (f^TTf) describes an unstable (upside-down) 
harmonic oscillator 

The classical solution for a mode in the system can be 
written as 



C V) = ^<i( fc )«i(^. 1) + 4>% S {k)u 2 {r,\ 77), (21) 



where 



ui = 



U-2 



sin[fc(?7 - + fcW) + cos[fc(ry - r)')](r)' - rj) 

sin[fc(r7' - rn)\{\ + krurf) + cos[fc(ry' - rji)](rn - 77') ' 
sin[fcfa - j)]{\ + krfrjj) + cos[k(r]' - rjj)](r]' - rg) 

sm[k(i]i - rj)](± + krjrji) + cos[k(r/ - rn)](r) - rji) 



We will assume that the system-field contains only one 
Fourier mode with k = fco- This is a sort of "minisu- 
perspace" approximation for the system-field that will 
greatly simplify the calculations, therefore we assume 

<l>^ cl (x,r,') = ^\r,')co S (k .x), (23) 

where 0i cl is given by 1211) . 

In order to obtain the master equation we must com- 
pute the final time derivative of the propagator J r . After 
that, all the dependence on the initial field configura- 
tions 4>%i (coming from the classical solutions </>< cl ) must 
be eliminated. Followingthe same procedure outlined 
in previous publications |9j, we can prove that the free 
propagator satisfies 



2a 2 (77)^(77, 77)^1(77' ,rj)<t>% } 



(24) 



a 2 (vWi(v,v) ~ a 2 (rn)u' 2 (rji,ri) 

M 1x1(77', rj)V ^ 
[a 2 (77)74(77, 77) - a 2 (r)i)u' 2 (r]i,r))] 



where a prime now stands for a derivative with respect 
to ?/ and the spatial volume V appears because of we 
are considering only one Fourier mode for the system. 
These identities allow us to remove the initial field con- 
figurations 4>%i, by expressing them in terms of the final 

amplitudes 6^ f and the derivatives d,± , and obtain the 

<-j <p <f 

master equation. 

The full equation is very complicated and, as for 
quantum Brownian motion, it depends on the system- 
environment coupling. In what follows we will compute 
the diffusion coefficients for the different couplings de- 
scribed in the previous section. As we are solely inter- 
ested in decoherence, it is sufficient to calculate the cor- 
rection to the usual unitary evolution coming from the 
imaginary part of the influence action. The result reads 

i 9.7Pr[^</|0</;»/] = (0</|[-ffren,Pr]|0</) 

- i r 1 D 1 (k ,ij,A)+r 2 D 2 (k ,ri,A) Pr#</| 

+ (25) 
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where we have defined T\ 



T2 = -^p(0<y — ) 2 . The ellipsis denotes other 
terms coming from the time derivative that not con- 
tribute to the diffusive effects. This equation contains 
time-dependent diffusion coefficients Dj . Up to one loop, 
only D\ and D2 survive. Coefficient Di is related to the 
interaction term </A 0> , while D2 to <j)\ (/)> . These coeffi- 
cients can be (formally) written as 



and fined as the solution to 



£>!$,, »?, A) = f \ 'a\T,)a\r,')F*(T,,T,' \k ) (26) 



and 



x ImGlAr], r)', 3fc ) 9(3fc - A), 



D 2 (k , V ,A) = -36 / drfa\ri)a\if)F^{ri,rf,ko) (27) 

Jrii 

x [ReGf + (r,, rf, 2k ) + 2 ReGf + (r,, rf, 0)], 



(30) 



fVd 

3 Jr ' i 
rVd 

> T ; / driDi(ka,A,v), 

Jr)i 



where the inequality is valid for any particular j = I. 
That is, the interactions with the environment have a 
cumulative effect on the onset of classical behaviour, i.e. 
the inclusion of a further interaction term reduces the 
decoherence time rjd- Therefore, in order to find upper 
bounds to 77^, we define the decoherence time 77^ . coming 
from each diffusion coefficient by 



1 w Tj / di] Dj(ko,A.,<i]), 



(31) 



with j = 1,2. 



with the function F c \ defined by 

j? 1 j \ sin[fc (T?- %)] Th cos[k (r] - r] t )} 
F c i(r],T]i, fc ) = — 1 . (28) 



koV 



The explicit expressions of these coefficients are compli- 
cated functions of conformal time, the particular mode 
fco, and the cut-off A, and we show them in Appendix lAl 
It is important to note that here we are only studying 
the effect of normal diffusion terms, even it is known 
that anomalous diffusion terms can also be relevant at 
zero temperature. Analysis done in Ref. (2^ suggests 
that anomalous diffusion for a supraohmnic environment 
is only relevant on a small transient and decoherence for 
unstable long-wavelength modes are driven by normal 
diffusion coefficients l28l. 



A. Diffusion terms: Numerical results and analytic 
appr oximat ions 

In this subsection we will analyze the behaviour of each 
diffusion coefficient as a function of the Fourier mode 
fco (considered for the system in Ea. H23|l ) and the cut- 
off A. We will also analyze the temporal evolution of 
the coefficients and their integration in time for fixed 
values of fco and A. We will present simple analytical 
approximations to the coefficients which can be used in 
Ea. (l31|) instead of the full expressions to estimate the 
decoherence time-scale. 

We define the dimensionless quantity Af [fco, rj\ which is 
the number of e-foldings between the time 7]k when the 
mode fco crosses the Hubble radius (i.e., |fco^fc | — 1) and 
any time rj during inflation, 



IV. DECOHERENCE 

Coherences are destroyed by diffusion terms. This pro- 
cess is evident after considering the following appoximate 
solution to the master equation 



Pr [^<, 0<; T7] ~ Pr 



dr\ £>j(fc , A, 77) 



(29) 



where is the solution of the unitary part of the master 
equation (i.e. without environment), and Tj includes the 
coefficients in front each diffusion term in Ea. H25|) . The 
system will decohere when the non-diagonal elements of 
the reduced density matrix are much smaller than the 
diagonal ones. 

The decoherence time-scale sets the time after which 
we have a classical field configuration, and it can be de- 



Af[ko,rj\ = In 



Vk 
V 



In I fc 7? I 



(32) 



This quantity has the special feature that its sign indi- 
cates whether the mode is inside (Af[ko, 77] < 0) or outside 
(W[fco,?7] > 0) the Hubble radius. 

Let us first consider the diffusion coefficient Dt, which 
comes from the interaction term 0^0>. Because of we 
are considering only one Fourier mode for the system, 
with wave vector fco, and the environment-field contains 
only modes with fc > A, this coefficient is different from 
zero only if A/3 < fco < A (i.e., <j> < is only coupled with 
the fc = 3fco mode of the environment). 

For this coefficient we can obtain an exact analytical 
expression from Eq. l|A3|) (see Appendix lA"|l . In Fig. Q] 
we have plotted this expression as a function of fco for 
a particular value of the conformal time (if 77 = —1/2). 
For later times, the graphs are qualitatively similar but 
D\ oscillates more rapidly (since we have obtained our 
results by perturbative calculations, they are not valid at 
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FIG. 1: Coefficient D\ (in logarithmic scale) as a function of 
k for H V = -0.5. 
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FIG. 2: Coefficient D\ (in logarithmic scale) as a function of 
fcn for a fixed value of |ftoJj| = 1/10 (feo > H/10), where we 
have also plotted the approximation in Eg. 13311 . 

large times). The coefficient decreases with fco and takes 
its maximum value for fco = A/3, implying that this kind 
of interaction produce more decoherence for small values 
of fco and, from the above discussion, for small values 
of A. For practical purposes, we consider the following 
simple approximation to D±: 

Dr r ° x (ko,V,A) = -Woll^M®^ A) - (33) 

As we can see from Fig. [21 this approximation is less 
close to the exact coefficient for big values of fco- It is 
important to note from Eq. (|3 l|t that if the approximation 
is a lower bound to the coefficient, then it will be useful 
to calculate an upper bound to 77^ . 
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N[k ,Tl d ] 

FIG. 3: Temporal integral of D\ in units of H~ 1 as a function 
of Af[ko,rjd] for two different values of ko- Solid curves are 
numerical integrations of the exact coefficient, while short- 
dashed curves are from analytical integration of Eq. 1331 . Note 
the logarithmic scale on the vertical axis. 

According to the definition of 77^, given fc , A and r 1; 
we can estimate this time by integrating D\ over the 
conformal time 77 and plotting this temporal integral as 
a function of A/"[feo>%i]- Fig- El shows such plots for 
two particular values of ko/H, where we have added the 
curves of the approximation (|33(l . As we have illustrated 
with these examples, the approximation to D\ is useful 
to estimate the order of magnitude of the corresponding 
decoherence time 77^ ■ 

Let us now examine the behaviour of the coefficient 
D2 which is associated with the interaction term (fy^c/y^. 
Since the interaction is now quadratic in (f> > , there are no 
restrictions on the values of fco such that D2 7^ 0. Hence 
this coefficient can affect the coherence of all modes in the 
system, therefore it is the most important in our model. 

The dependence of D 2 with A is showed in Fig. 0] 
for three different values of fco and a fixed time 77 = 
— 10~ 3 H~ 1 . As we can see in this figure, for modes with 
fco < < A the coefficient is weakly dependent of the value 
of the critical wave vector A. 

In Figs. |S] and E] we have plotted the coefficient D2 as 
a function of kg/H for fixed values of 77 and A. Fig. El 
shows that modes within the Hubble radius (\kor)\ > 1) 
the diffusion coefficient is an oscillatory function and it 
has a maximum when fco ~ A. The same behaviour was 
noted for conformally coupled fields [24| . Physically, this 
fact can be interpreted in terms of particle creation in 
the environment due to its interaction with the system. 
For these modes, a very simple but good approximation 
to D2 is given by 

07 1 

where the upper-script £ stands for "local" , since it corre- 
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FIG. 4: Coefficient D2 (in logarithmic scale) as a function 
of A/H for a fixed conformal time r) = — 10 _s i? _1 and three 
different values of fc (H/W, H and W0H). 




FIG. 5: Coefficient D2 as a function of fco for A = 20H and 
Hr) = -1/2. 



sponds to approximate the whole diffusion coefficient by 
the term in Ea. 1)27(1 containing the Dirac delta function 
(see details in Appendix 0). On the other hand, Fig. 
shows the dependence of D 2 with fc for modes outside 
the Hubble radius at 77. For these modes we note that 
the diffusive effects are more important for the smallest 
values of fco, which are most sensitive at the expansion 
of the universe. Note that the classical equation 12U|) 
for each mode of the free field ip (defined as ip — a(r])(f>) 
describes an stable (unstable) oscillator if the mode fco 
satisfies |fco?7| >> 1 (\koT]\ << 1). Therefore, it was ex- 
pectable that the diffusion coefficients mirrors this fact. 

For modes far outside (\k ri\ « 1) the Hubble radius 
we have found an asymptotic approximation of D 2 , useful 



FIG. 6: The same as Fig. |K|but for modes fco outside the 
Hubble radius at rj = — l/2jy _1 and A = 2H. 




FIG. 7: Coefficient D2 (in logarithmic scale) as a function 
of Af[k ,rj d ] for A = 300H and fc = 100#. Heavy solid 
curve is numerically calculated from Ea. l|27^ (see details in 
Appendix Light solid curve corresponds to the approx- 
imation obtained by averaging the local (dotted curve) and 
the asymptotic (short-dashed curve) approximations. 



if A ~ fco j which can be written as 



1 



7T 2 iJ 4 AV ' 



(35) 



Fig. [7| shows a graph of the coefficient D 2 as a func- 
tion of A/"[fco, ?7ci 2 ] for a particular value of fco and A, and 
the curves of the two approximations above. From this 
figure, we can distinguish two different regimes: one is 
well described by the local approximation (|fco?7| > 1) 
and the other by the asymptotic one (|fco?7| << 1). This 
behaviour of D 2 indicates that the decoherence process 
becomes faster once the mode crosses the Hubble radius. 
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FIG. 8: D2 and D^pprox ag f unc ^j ons of fc /// for a fixed value 
of |fc r?| = 1/10 (fc > H/10) and A = 100//. The vertical axis 
is on a logarithmic scale. 
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FIG. 9: Temporal integrals of D2 and ]j^PP rox ( on a logarith- 
mic scale and in units of // _1 ) as functions of Af[ko,r/d], for 
A = 300// and two different values of k (100// and 200//). 



In addition, we can see that a better approximation to 
D2 is obtained by averaging the local and the asymptotic 
ones, that is 



j-.approx ^2 



(36) 



which is also showed in the same figure. 

The dependence of D2 with ko for a fixed value of 
Af[ko,rj\ is shown in Fig. [H] in which we see that the 
approximation is less close to the numerical curve for big 
values of ko, but it bounds D2 from below. 

In Figs. IH1 and ^| we have plotted the temporal in- 
tegral of the coefficient D2 as a function of Af[ko,r]d 2 ] 
for different values of ko and A, where we have also plot- 



FIG. 10: The same as Fig. [§] but for smaller values of ko 
(k = H and k = H/10) and A = 2H. 



ted the curves obtained using the approximate expression 
given in Ea. l|3tj|) . In view of our analysis, it is reasonable 
to use that approximation to estimate the decoherence 
time-scale r\d 2 , at least for values of A smaller than 10 3 fco. 



B. The decoherence time 

So far in this Section we have analyzed the behaviour 
of the diffusion coefficients, and we have shown the ap- 
proximations considered are useful to estimate the de- 
coherence time associated with each diffusion coefficient 
(or interaction term) . Let us now work at the level of the 
order of magnitude and apply these results to quantify 
the decoherence time rid- 
As we have already mentioned, r)^ and r\d 2 are upper 
bounds to rjd- The time-scale given by rjd(ko) sets the 
time after which we are able to distinguish between two 
different amplitudes of the Fourier mode ko within the 
volume V . Thus, the maximum value of rjd(ko) with ko < 
A corresponds to an upper bound to the decoherence time 
for the system-field (since, in principle, <f> < contains all 
these modes with amplitudes different from zero). 

In order to quantify decoherence times 77^ 2 we have to 
fix the values of Ti^ (i-e., we have to assume values to A, 
V, and 0<A For this, as a first approximation and 
since we are considering a fixed de Sitter background, we 
will assume that the slow-roll conditions are satisfied at 
least up to times of the order of the decoherence time- 
scale. We will choose typical values for the parameters 
of the model (A, V) and for the elements of the reduced 
density matrix (0^,0~y). 

The slow- roll conditions are usually written as [29ll30| : 



16tt 



U" 
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where m p i = G -1 / 2 is the Plank mass, U — A</>g, 
V = dU/d<j)o, and 6 is the classical inflaton field. If we 
also assume that So is homogeneous on physical scales 
A p h >> (VA^o) 1 ! these conditions imply that the clas- 
sical configuration field do and the Hubble rate satisfy 



H 



&ttU 
3^ 



d<t>o 
dt 



EL 

3H' 



(38) 



Defining the end of the inflationary period setting 
eu ~ 1) one can set ^(N^) = W (N v + l)/ir m p i, where 
N v = lna(m)/a(?7); typically the e-fold number N Vi 
X fit) |29j| . Thus, we assume the mean value of the 
system- field at time of decoherence is 0o(-^T) d ) C^<p <f = 
(</»</ + 0</)/2 — 0o(^j) d ))- Using this mean value and 
Eq.(|2EJl we can write H 2 ~ &\N 2 d rri^ x / (?>-k) . In order to 
choose a typical order of magnitude for the fluctuations 
of the system-field, let us use the fact that the amplitude 
of the so-called primordial density perturbations 5 can 
be inferred to be of order 10~ 5 , and that 6 ~ H/doSd, 
where 56 is the amplitude of the inflaton field fluctua- 
tions [23, OH- Thus, with the use of these constraints and 
Eq.® we set = (6+ f -cj>- f ) ~ lQ-^6 (N Vi )/N Vd . 
Since V is the spatial volume inside which there are no 
coherent superpositions of macroscopically distinguish- 
able states for the system-field, it is reasonable to choose 
V = vH~ 3 , with v ~ 1 (which corresponds, if N ~ 60, to 
the comoving volume V ~ (ao-ffo)~ 3 , where Hq and do 
are the present values of the Hubble rate and the scalar 
factor, respectively) . 

From previous considerations and according to the con- 
dition {SH , the times % 12 can be estimated as the solu- 
tion to 



where t c i 1 = H^ 1 hxa(rid^), and we have defined ko = olH. 

With the use of the approximation in Ea. l|36[) we obtain 
a quadratic equation in x = a 3 (rjd 2 ) = exp(3Htd 2 ), which 
is simple to solve for td 2 ■ The result is 




where a = A/H. For the sake of simplicity let us set 
a ~ 1. Since H/T 2 » 1, it yields 



d-2 



< J_ lnfl27rV — 
~ 6H \ r 2 



~ i? + 6i? n I v V 60 , 



(42) 



Assuming A = Ht cn< \ > 60 as an estimative scale to the 
end of inflationary period and values of A < 10~ 5 , we 
obtain 



iend ~ 60 



120 



hi 



(43) 



which makes sense only if ko < A < 3/co (a < a < 3a), 
and 

i ond 15 120 W 3 ) v ; 

From scales td x and i<j 2 we conclude that if one set 
A < H, the decoherence time-scale for the system-field is 
shorter than the minimal duration of inflation for all the 
wave- vectors within the system sector. 



^ l( ,)^^<2xlO- (A1 ° 5) ^ 



H 

r7 



v V 60 



— — (39a) 



HdrjD^) ~-<9xl0 17 - (-) , (39b) 
1 2 v \ 60 ' 



where we have used the fact that N Vd < N. The last term 



A *</ S *</- 



in (|39a|l follows after taking (6^ 3 - 6^) 
For example, from Fig. ^|we can see that the temporal 
integral of D2 takes a value of order 10 17 for Af[ko = 

H, r)d 2 ] between 7 and 8. As D2 is weakly dependent 
with A for ko < H, we can say that it is valid for A ~ H . 
Thus, since D2 decrease with kg for \koT)\ < lj we can 
conclude that the decoherence time for those modes with 
kg < H are smaller. 

Substituting the approximation given in Eq.© into 
the left-hand side of l|39a[l and assuming that l-ffjjdj << 

I, we get 



V. EFFECTIVE DYNAMICAL EVOLUTION OF 
THE SYSTEM-FIELD IN A FIXED DE SITTER 
BACKGROUND 

In this Section we concern ourselves with the time 
evolution of the system-field. After reviewing the phe- 
nomenological way to describe the stochastic dynami- 
cal evolution of the system, we present the renormalized 
"semiclassical-Langevin" equation. This equation can be 
used to describe the dynamical evolution of the classi- 
cal configurations of the system-field and hence is useful 
once all the modes in the system have lost coherence. As 
a first step to understand the generation of classical inho- 
mogeneities from quantum fluctuations, we then consider 
a simple situation in which we analyze the influence of 
the environment on the power spectrum for modes inside 
the system sector. 



< 
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A. Effective dynamical evolution, noise and 

(40) 

expectation values 

The real part of the influence action contains divergent 
terms and should be renormalized. The imaginary part is 
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finite and is associated with the decoherence process. It 
is well known that the terms of the imaginary part that 
come from a given interaction term in the original action 
can be viewed as arising from a noise source 0, . In 
our case there are two such sources £2 and £3, which are 
associated with the interaction terms and </><</>> 

respectively. That is, the imaginary part of the influence 
action can be rewritten as 



ImSA=-\n(F[A 2 ]F[A 3 ]), 



(45) 



where i^[A„] (n = 2, 3) is the characteristic functional of 
the noise £„, which is related with Gaussian functional 
probability distribution -P[£ n ] as 

F[A n ] = J D£ n P[£„]expj-i j d A xA n {x)^ n {x)^ , 

p [$n] = N„exp{-i J d^Xi Jd 4: X2^n(xi) 

x is- 1 (x 1 ,X2)Z„(x 2 )y (46) 

where N n is a normalization factor and v~ l is the func- 
tional inverse of the noise kernel v n : 

v 2 (x u x 2 ) = -288\ 2 a A { m )a A { m )ReGf + {x 1 ,x 2 ) 1 (47a) 
v 3 (xi,x 2 ) = M\ 2 a\ m )a\ m )ImG^ + {x u x 2 ). (47b) 

The Gaussian noise field S, n (x) is completely character- 
ized by 

(£n(x)) P = 0, (48a) 

{€n(?l)€n(X2))p = Vn(xi,X 2 ), (48b) 

where with ()p we are denoting average over all realiza- 
tions of £ n (x). 

The functional variation 



SS C 



= o, 



(49) 



yields the "semiclassical-Langevin" equation for the 
system-field, which is only valid once the system has be- 
come classical. 

With the identifications above, the reduced density 
matrix can be rewritten as [32^: 

pA<t>% f \<t><fM = j DfrPfa] J D&P[&] (50) 

with 



Pr[0< / |0 </ ,6,6;?7] 



X / D(j) < Pr[^><i|0 <i ; Vi] ex P {iSefi[4»<, £2, £3]} 



where the effective action S e s is given by 
S eS [<P+, = Re{A[^+A<]} (52) 

d 4 x[A 2 (x)&{x) + A 3 (x)&(x)}. 

Here A is the CGEA specified in Ea.fnty. Thus, the full 
expectation value of any operator Q[<^><] can be written 

as 



(W<]> = / D&Pfa] j D&P[&] J <#< 

x Pr[(t><\<f><,&,&,rj\ Q[(t><] 

= ({Q[<t><]) q }p, (53) 

where Q q is an usual quantum average for a system-field 
subjected to external stochastic forces. 

B. Renormalized equation of motion for the system 
Taking the functional variation as in Eq. 149(1 we obtain 

<t>< (v,x)~ A0< fax) + 2 H<l> < (77, x) 

+ 4 Aa 2 {rj)4>< (77, x) [4> 2 K fax)- 3iG^ + (77, 77, 0)] 

-96X 2 a 2 fa) ( dj/aV)^^) 



288 X 2 a z fa / drj'^fa)^^^) 



x / d 3 ycj) 2 < (r]\y)ImGX+(v,v' ,x - y) 



A2 /„ / 



-&(ri,x 



a 2 fa 



3, , ^\fax) 



(54) 



This equation contains divergences. In order to renormal- 
izc it we use the method of adiabatic substraction with 
dimensional regularization, which works at the level of 
the field equation and is particularly useful for solving 
the equation numerically |33j. To simplify the task we 
assume the system-field to be homogeneous enough so 
that the spatial derivatives and the term in Eq. (|54|l com- 
ing from the interaction 0></>< are negligible. Note that 
for an homogeneous system-field there is no contribu- 
tion from this interaction term up to one loop order, due 
to orthogonality of the Fourier modes (see, for instance, 
Eq.JTJJ). Details of the renormalization procedure are 
relegated to Appendix FBI The renormalized equation is 



^(77) + [AM 2 fa + AZfaR] a 2 fa Mv) 
+ 2 H ^fa + 4 [A + AXfa] a 2 (77) ^(77) 

+ -—^faa'fa / drf 4> 2 <{v')Jfav') 



(55) 



" , , 4><WW<W) „, , u < 

11 Halrf) ' = 



a 2 (77) 
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where R = 12H 2 ; we have redefined the noise source as 
1 



6 fa) = y J d 3 x&(v,Z), 

whose correlation function is given by 

36A 2 



(56) 



3[2A(7 ?1 - 7ft)] 

A 



1 



-(7r-2<Si[2A|ryi- 77 2 |]) 
- sin[2A(?7i - 772 



Vim 3 rftri 3 A 2 r7?r?| 



3 ?7 2 ?7| 



+ 



/71772 



1 



2(r?i - 



.2(th-i») 3AV% 2 
and we have also defined the following functions: 



9A 2 
2tt z 



1 


-In 


a{ri)n 




2 


2A 




"7 


-In 


a(r])fi 




2 


2A 





AM^^A^-^^fe) 

x Ct[2A(n - tk)] 
V 



27r 2 a 2 fa)' 

H{ V ,rf) = -^Ci[2A{r,-rf)], 
77' 

J(V, rf) = Ml, ff) + J2 (V, rf) + Mv, rf), 
J 1 (r 1 ,r l ')=ei[2A(i 1 -r l ')} 



2fa 3 -»n , r? 2 

/3 



3r/ 4 

3?y' A- 
sin[2A(r/ - 77')] 



'/ 



^3(r?,V) 



2A(77 - 77') 
sin[2A(?7 - 77')] 
A 



3;/ 



,/4 



+ 



(57a) 
(57b) 
(57c) 

(57d) 
(57e) 

(57f) 
(57g) 

(57h) 



where Ci[x] is the cosine integral function |3J]. Notice 
that these functions are logarithmically divergent in the 
limit A —> 0, which is the well-known infrared divergence 

M 

The functions above contain useful information. Par- 
ticularly, they allows us to examine the conditions under 
which the loop expansion breaks down. In order to esti- 
mate the time after which the one-loop terms become of 
the same order of magnitude as the classical ones we may 
compare their time-dependent parts. For example, from 
Ea. (|57ll we see that the time-dependent part of AA(r/) 
(AE(r?)) is the order A (one) for Ht~l/\ and AM 2 (77) 
is important only at the initial time (77 ~ rji). In the same 
way, we can use the term Aa 2 0< to compare it with the 
remaining ones. In Fig. II II we have plotted the Ji func- 
tions, where we can see that J\ dominates all others. It 
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FIG. 11: J and Ji. 2 , 3 as functions of Af[A, rj'] = - In |Ar/| for 
At? = -KT 2 (AT [A, rj] = In 100 ~ 4.605). 



allow us to make the following approximation: 



1 



-^(^HAt,!) 2 , 



(58) 



where the last term is a simple long-time asymptotic ex- 
pression. With the use of this approximation we obtain 
that the two terms in question are of the same order when 
|ln|A?7|| - ITS/A and, for A ~ H, Ht ~ 1/VA. 
As it is shown in Fig. ^1 the TZ function peaks at 

77 - 



77' and hence we can approximate 

. Ha(rj) X Ha(ri) J m 

Ha{rj) 
x ( 7 + ln|2A77|), 



(60) 



where the last expression corresponds to a long-time ap- 
proximation. As for the J term, using the long-time 
approximation we get In \Arj\ ~ (j) < H/(\(j) < ). If in addi- 
tion we assume that the slow-roll conditions are satisfied, 
this time can be estimated as | In 1 7V77 1 1 ~ l/(Xeu) and, 
for A ~ H, Ht ~ l/(\eu), where ejj is the slow-roll 
parameter given in Ea. l|37|l . 

From previous discussion, we can see that the addi- 
tional terms are only important at times longer than the 
typical time-scale associated with the decoherence pro- 
cess. 



C. Generation of inhomogeneities: role of the noise 

Having reached this point, we ask ourself about the 
primordial inhomogeneities. Certainly, the model we are 
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FIG. 12: Kasa function of N[A,r)'] for At? 
{N[A., n] = In 100 ~ 4.605). 
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considering is too simplified to obtain any property of 
the power spectrum of such inhomogeneities. Neverthe- 
less, we can discuss some relevant aspects about the way 
of connecting the initial quantum fluctuations with those 
which eventually become classical and statistical. With 
this purpose, let us compute the power spectrum for some 
of the Fourier modes in the system working in an approx- 
imate way and making certain assumptions, but taking 
into account the decoherence process. 

According to the results of the previous Section, the de- 
coherence process occurs at different rates for each mode 
in the system. Provided that decoherence is more effec- 
tive for modes outside the Hubble radius, we consider 
the following situation. We suppose that the "relevant" 
modes are in the Bunch-Davies state (i.e., in equilibrium 
with the environment) at the initial time r\i. Conse- 
quently, the matrix elements of the initial density op- 
erator for each of these modes are Gaussian functions in 
a mode- amplitude basis. As a first approximation we ne- 
glect the non-linearities which might affect the Gaussian 
form of those matrix elements. 

We are now interested only in computing the power 
spectrum of those "relevant" modes up to h and A 2 
order. To carry this out, we split the system-field as 
4> < = 4>a{ii) + S4> < , where we identify (froiv) as a classi- 
cal background field which satisfies slow-roll conditions. 
The power spectrum of the quantum field fluctuations 
5(f> < may be expressed as P^(k) = 2ir 2 k~ 3 A^(fc) , with 
A 2 (fc) defined by 



: (x)<^<(£+r)) = / d 3 k 



A|(fc) 
Ank 3 



exp(— ik ■ r). (61) 



where () is the full average in Eq. I|53|1 . 

Consistently with the assumption that the "relevant" 
matrix elements remain Gaussian, we expand the semi- 
classical equation (|54|) up to linear order in the mode 



amplitude of interest 8<fi<(k). 
we obtain 



Through this procedure 



^'(7 ? ) + 27^(7?) + 4Aa 2 ^( r? ) =0 



8^{k, V ) + [k 2 + 
+ 2H8cj>' < (k,r 1 ) = 



12AaVo«<MM) 
2 — VoW), 



(62a) 



(62b) 



where we have discarded the terms which do not con- 
tribute to the power spectrum up to h order. The term 
with the £ 3 noise source gives a zero contribution due to 
our approximations and the orthogonality of the Fourier 
modes. It is important to note the presence of the £2 
noise source, which ensure the decoherence process oc- 
curs as the matrix elements are evolved in time. 

In order to obtain the power spectrum let us split 
(5</><(fc, 77) = <5</><(fc, 77) + Stfil- (k, 77), with 50^-(fc, 77) = 
(5cj> < (k,ri))q, where () 9 is the quantum average defined 
in Ea. (|53|l . Because of the assumption of linearity, 
this quantum average satisfies the semiclassical equation 
(|62b(l . whose solution can be written as the sum of the 
homogeneous solution 8<f> 1 ^(k, 77) and a particular solution 
<5</><(fc, 77). The former is given by 

8<fe(k,ri) =a- 1 (77)Kv^l^ + /3fc^|J^], (63) 
where otk and [3k are constants of integration, and v = 
6\<fi 2 /H 2 . Setting v ~ |, a particular 



. — e, with e 
solution is 

where 



dv'9(k,V,v')&(k,v')?o(v'), (64) 



sin[/s(77 — 77')] 



0(77)0(77') 
cos[fc(77 - 



k 2 f]rj' 



'■r\r\ 



(77 - 77') 



(65) 



With the use of the initial conditions (o^<(fc, rji)) q = 
(<5</>< (k, rji)) q — 0, we obtain that (<!></>< (fc, 77)) p = 0, and 
thus 8(f>i : (k,ri) = (£;, 77). Within these approxima- 
tions, the result is analogous to that for the linear quan- 
tum Brownian motion (QBM) [2?], HU. Therefore, the 
quantity A|(fc) can be written so that it receives two 
contributions: 

A 2 (k) = A 2 q (k) + Al 6 (k). (66) 

The first one comes from the unitary evolution of the 
initial density matrix, i.e., it is the usual quantum result 
for the case of the free field ;29]: A 2 ,(fc) = (H/2n) 2 (1 + 
k 2 rf). The second one appears due to the £2 noise source 
and can be computed through 

{84> 6 < {ki,ri)8(j)^{k2,ri))p = f dru f d^oiv^Voim) 

J T)i J Tji 

x g(ki,v,vi)g( k 2,v,V2), (67) 
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FIG. 13: A^/(A(£o) 2 as a function of Af[k , V ] for A = 5H 
and two different values of ko- Note that the vertical axis is 
on a logarithmic scale. 



where 

(6(fei,??i)^(fe,r 72 ))p = -288A 2 a 4 (r 7 i)a 4 (r? 1 ) 
x(27r) 3 <5 3 (fc; - k 2 )ReGf + {n um M). (68) 

Thus A 2 ^ (k) can be expressed as 

144 F> r n 
A^ 5 (fc) = -A 2 — fc 3 / / ^ 2 a 4 ( % )a 4 ( % ) 

x (j> (rn)(t)Q(v2)g(k., r), m)g( k , v, m) 

x ReGf + { m , m ,k). (69) 

The Fourier transform R eG^ _(r]i, 772 ; fc) (where fc < A) 
can be obtained from Ea. (|A9(l (see Appendix replac- 
ing 2fc by fc. 

Since the additional contribution A 2 ^ (fc) to the power 

spectrum is of order A 2 , it is expected to be negligi- 
ble. As we are assuming that 4>q is a slowly varying 
field, we can compare the relative order of magnitude 
of both contribution setting <f>o(r)i) ~ </>o( ? ?2) ~ <fio an d 
H 2 ~ A^o/m 2 ,}, so that <fo can be taken out from the inte- 
gration in Eq. (|69|l . Thereby, when \krj\ < 1, the contribu- 
tions A^ 5 (fc) is negligible compared to the usual A 2 , (fc) 

if (A^£ (fc)/A0 o ) 2 << A -1 (0o/m p i) 2 . Since the last quo- 
tient is usually much bigger than one, this condition is 
typically satisfied if (A^; (fc)/A</>p) 2 < 1. As in the ex- 
ample shown in Fig. this is the case and hence the 
additional contribution can be neglected. 

On the other hand, the usual contribution A|, (fc) is 
independent of fc for a fixed value of kr], corresponding 
to a nearly scale-invariant spectrum, whereas A 2 ^ (fc) de- 
pends on fc and A (see Figs. H4land lT5|l . To see this more 
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FIG. 14: A^/(A0o) 2 as a function of k /H for A = 5H and 
\k v\ = 0.001. 

0.16 1 1 1 1 1 1 1 1 1 1 1 




2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 
A/H 



FIG. 15: A^/(A<£ ) a as a function of A/H for k Q = 2H and 
|fco?7| = 0.001. 

clearly, it is useful to rewrite A 2 5 (fc) as 

A2 m 144 [ kVdx ^ [ kVdx ^ ( U < n \ 

A ^( fc ) = / -zr / -rMvMm) 

fc 

where / = k 3 Hg and F = k z H~ l ReG t ^ 2 + , which de- 
pend on fc not only through the kr\ combination, but also 
through the combinations krji and A/fc. It is well-known 
that a finite duration of the inflation stage produces a de- 
parture of the power spectrum from the scale-invariant 
one . The breaking of the scale invariance by the pres- 
ence of the infrared cut-off A is similar to the one found 
in Ref. [H. 

It is important to note that we have assumed a Gaus- 
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sian initial condition for the reduced density matrix ele- 
ments and that the Gaussian form of these elements is not 
affected by the non-linearities of the interactions. How- 
ever, as was pointed out in Ref. ^(Jj ^ is expectable that 
when the non-linearities become important, the Gaussian 
approximation breaks down and therefore the associated 
Wigner functional becomes non-positive. Nevertheless, 
as well as for a given mode the Gaussian approximation 
remains valid up to times longer than the decoherence 
time-scale, one can use the classical description for it 
even in the non-linear regime. 



tion plus the one which appears because of the system- 
environment interaction, we found that the latter is neg- 
ligible compared to the former. In spite of this result, we 
have remarked that the system-environment interaction 
is essential to have a complete quantum to classical tran- 
sition which allows a late-time classical treatment of the 
system degrees of freedom. 

As for future work, we consider that it is worth apply- 
ing the same procedure used in this paper to inflationary 
models involving two or more interacting fields, such as 
some hybrid model [3^| . 



VI. FINAL REMARKS 



VII. ACKNOWLEDGMENTS 



Let us summarize the results contained in this paper. 
After the integration of the high frequency modes in Sec- 
tion II, we obtained the CGEA for the modes whose wave 
vector is shorter than a critical value A. From the imag- 
inary part of the CGEA we obtained, in Section III, the 
diffusion coefficients of the master equation. System and 
environment are two sectors of a single scalar field, and 
the results depend on the "size" of these sectors, which 
is fixed by the cut-off A. 

To analyze the decoherence process, in Section IV we 
evaluated the diffusion coefficients and its integrations 
over the conformal time. It allows us to conclude that 
if we consider a cut-off A ~ H , those modes with wave 
vector fc -C A are the more affected by diffusion through 
one of the coefficients. For these modes, we saw that the 
effect is weakly dependent of the critical wave vector A. 
If one consider a cut-off A > H, and modes H < ko < A, 
diffusive effects are larger for those modes in the system 
whose wave vector is close to the critical A [24| . 

We presented the complete expression of the diffusive 
terms (in Appendix A), and also some simple analytical 
approximations to the coefficients, which are useful to 
make an evaluation of the decoherence time-scale. We 
performed an extensive analysis of the evaluation of the 
time-scale for the decoherence process for a typical case 
of interest. In such a case, we obtained that for a given 
mode fco < A < H, decoherence is effective by the time 
in which inflation is ending. 

In Section V we analyzed the effective evolution of the 
system. Assuming an homogeneous system-field we first 
presented the explicit form of the renormalized stochastic 
Langevin equation and we analyzed the relative impor- 
tance of the terms appearing in that equation due to the 
system-environment interaction. From such analysis we 
conclude that those terms are of the same order of mag- 
nitude than the classical ones for times longer than the 
typical decoherence time-scale. 

We them considered inhomogeneity generation in a 
simple particular situation, in which we analyzed the in- 
fluence of the environment on the power spectrum for 
some modes in the system. In that situation, splitting 
the spectrum as the sum of the contribution coming from 
the unitary evolution of the Bunch-Davies initial condi- 
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APPENDIX A: DIFFUSION COEFFICIENTS 

In this appendix we describe some technical details 
about the computation of the diffusion coefficients Di 
and D 2 , starting from Eq.(J5SJ) and (|T7Jl respectively. 

In order to evaluate the coefficient D\ , we have to per- 
form the Fourier transform of the imaginary part of the 
propagator for the environment G+ + (see Eg. 11260 ). This 
propagator can be expressed in terms of the mode func- 
tions (f>k given by 

1 exp{-ikr]} ( i \ -. 
Vk (2^) 3 / 2 V2ka( V ) \ krjj ^ S 



(j>k(r)) exp{ifc • x}, 



(Al) 



which corresponds to the Bunch-Davies vacuum assumed 
for the environment-field. The imaginary part of the 
propagator reads 

ImGl + (x 1 ,x 2 )= [ d 3 kRe{4> n ( Xl )(l>t(x2)}. (A2) 

Jk>A 

Substituting this expression into Ea. (|26|l we obtain 



dy 



cos [3 A 



■ryl 



1 



1 

9xy 



with 



F c i{x,y) 



sm[3Aag]Agg 
lfeV 

sin[A xy ] 



6a; 3 y 3 

F*(x,y) 9(3 -A), (A3) 



y cos [A 



■ryl 



(A4) 



x x 

where to reduce the notation we have defined the follow- 
ing set of dimensionless variables: 



H 



A 

fc 



V = I Wl i x = \k Q rj\ , A xy =x-y. 
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The integrations above can be performed exactly and the 
result can be numerically integrated over the conformal 
time 77. 

To obtain D2 from Eq. l|27|l we need to compute the 
Fourier transform of the real part of in k = 2ko and 
k = 0. This real part is given by 

Re G+% (77,77', 2k ) = -(2tt) 3 f d 3 k [ d 3 k' 

Jk>A Jk'>A 

xS 3 {k + k'- 2fc ) Re{^ k (r])fyv')4>k'(vWk'(v')},(^) 

where 4>k is the time-dependent mode function defined 
in Ea. HAlfl and <j>^ its complex conjugate. This Fourier 
transform can be expressed in terms of integrals of the 
form, 



1 m+n— 3 

T c _ K o 



d 3 k 

k>A Jk'>A 



d 3 k'S 3 (k + k' -2k ) 

.1 •-, A \ ' 



cos[(fc + k')(i] - r)')] 





k n k' m 




r+oc 


du r+ 2 


dz 


J A+2 


U"- 1 Ju-2 


z m-l 


/■A+2 


du r+ 2 


dz 


J A 


U 71 ' 1 J A 


z m-l 



where n and m are integer numbers (only m — 3 result 
to be necessary). The second equality follows after the 
change of variables u — k^ 1 k and z = k^\k — 2k$\ = 
ky 1 \jk 2 — 4fcfco cos(#) + 4fcg, where 9 is the angle be- 
tween k and fco- We also define the integral as the 
one obtained from I^ m by replacing the cosine function 
by sine. 

An integration of I^'m by parts yields 



n ' m + xy 2 - m 



+ °° du cos[2(u -I- l)A xy ] 
A m"- 1 (2- m)(u + 2) m " 2 

+ °° du cos[2(u - l)A xy ] 
a+2 u"- 1 (2-m)(A + u)™- 2 

A+2 du cos[2(w + l)A I 1 



A u" 



-,(A7) 



P - A 



n,m — 1 

2 - m 



(2 - m)A r ' 



sin[2(u + ^A^,,] 
A (2-m)(u + 2) m ~ : 

+°° du sin[2(u- l)A X j / ] 
A+2 (2-m)(A + u)™- 2 

A+2 du sin[2(u + l)Aaj. 1 
a ~ n 



(2 - rn)A m - 2 



(A8) 



With the use of these properties and definitions, the 
Fourier transform l|A5(l becomes 



where we have defined: 



(AlOa) 
(AlOb) 

c =J 3 3 - ^^32 - ^^22 (AlOc) 
These last integrals are easily computed, with the result 



Ib —I31 A-xyl^l 

1(1 =I?1 — 2A»,,Iqo — A 2 Jo 



I A = 2,6 (A xy ) + C ° S[2(A + 1)A 7J ~ C0S[2AA ^ ] ,(A11) 



A 2 



I B = cos[2A^]{Ci[2(A + 2)|A xy |] -Ci[2A|A KJ/ |]} 

- sm[2|A XJ/ |] {7r-«Si[2(A + 2)|A X3/ |]-«Si[2A|A XJ/ |]} 
sin[2(A + l)A xy ] sin[2AA xy ] 



AA 



xy 



AA 



(A12) 



Ic = {7T-Si[2(A + 2)|A xa |]-Si[2A|A x?/ rj} 

x (cos[2A^]|A x ,|-^^) 

+ {ci[2(A + 2)|A xy |]-Ci[2A|A x?/ |]} 
cos[2A X!/ ] ' 



x (sin[2A a .,,]A a .j, - r — 

cos[2AAa. y ] cos[2(A+ l)A xy ] 



A 2 



A 



(A13) 



where Si[x] and Cz[x] are the sine and the cosine integral, 
respectively [34| . 

The Fourier transform in k = is much simpler to 
compute, and the result is 



sin[2AA xy ] ( 


' 1 


f 2 A xy \ 


2A 


3 A 2 x 2 y 2 J 


cos[2AA xy ] 


(" + 


A 2 1 


A 


\xy 


3x 2 u 2 3A 2 x 2 y 2 



- (tt - 25z[2A|A xy |])^ + X ) } ( A14 ) 

Substitution of Eas. (|X9)l and l|A14l) into Eq.lO yields 
an expression for D2 that allows to evaluate it and its 
integration over the conformal time 77 numerically. 



APPENDIX B: RENORMALIZATION 

In this Appendix we obtain the renormalized semiclas- 
sical equation for the system-field given in Ea. (|55|l . To 
do so, we follow the same procedure as for the renormal- 
ization of the evolution equation for the mean value ((/>) 

m. 
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Starting with the bare action for the field (j>, the semi- 
classical equation for the homogeneous system-field </>< 
reads 

<P<(v) + H + Z R}a 2 (r))<t><(v) + W<(»7) 
+4A a 2 (77)^(77) - m a 2 (77)^ (77)^, (77,77,0) 



288AV (77)^(77) / drfa\rf) / tfytpM) 



4><{v) 

a 2 (r]) ' 



XImG ++ (r},r)',x - y) = -&{r},£) 



(Bl) 



[m = 0, £ = and A), ft is important to notice that only 
after subtracting the bare parameters can be replaced 
by the renormalized ones, since one of the integrals over 
k < A has an infrared logarithmic divergence. 

Writing the bare parameters in terms of the renormal- 
ized ones plus conterterms, 



+ Sm 2 , £ = + 6£ , A = A + SX, 



we can see that the divergences appearing for n — > 4 are 
cancelled with the use of the following counterterms: 



where mo, £0 and Ao are the bare constants, and R = 
12H 2 . Here we have replaced Ao by the renormalized 
constant A in the one-loop terms. 

Let us add and substrate the following terms to the 
left-hand side of Ea. lfBl")) : 



12Aa 2 (77)<M^)(0>)ad2 



(B2) 



where (4> 2 > ) a d2 is the adiabatic expansion of the expecta- 



tion value (</>>) up to the second adiabatic order 33]. 

In what follows we will show how the infinities in 
Eq.jBdJ are cancelled by the terms subtracted. We will 
compute the divergent terms added via dimensional reg- 
ularization, which will thus be able to be absorbed in the 
bare parameters as usual. 

The expansion of up to the second adiabatic order 
yields: 



with 



t) >)ad2 



where 



(B3) 



(B4a) 



(B4b) 



+ a 2 m,Q, and the integrations over the 



(4>l)ad2 = (4>l 


) F ad2 + < 


4>l) 


ad2i 


a f d 3 k 


fj__ 


(& 


-\)a 2 R 


2 J k>A (27ra)3 ] 








^[a' 2 + aa"} 




[aa 








a f d 3 k 


6A a' 






2 J k>A (2^a) 









wave vector k are restricted by k > A. 

In order to use dimensional regularization, we perform 
the integrations above over all wave vectors with k > 
and then we subtract the ones restricted by k < A with 
the result 



5m 2 =Am 2 , 



<5£ 

sx 



X 

: 4^ 



9A^ 
2tt 2 



n - 4 



AA, 



(B6a) 
(B6b) 

(B6c) 



where Am 2 , A£ and AA remain finite as n — > 4. 

Replacing the bare parameters by the renormalized 
ones in the integrals over k > A of Eq.|B4j) we obtain 



b>)ad2 



^>)ad2 



-iG* + (r),r),0), 



3A 
2^ 



fe>A 



dk 

T 



(B7a) 
(B7b) 



From equation above it is simple to note that the order 
A contribution in Ea. ijBlfl is completely cancelled by the 



one with (</>?, )f d2 (see Eq.{HU) 



In order to separate the divergent part from the order 
A 2 contribution of Ea. l|BlJl . we write the propagator as 



d*yImGf + (r,, rf, x, y) = ~(2tt) 3 / d 3 k Im[4>l 4>f 



dk ( cos[2fc(?7 — 7/)] 



fe>A 



a 2 (r])a 2 (r/^ 



fe>A 

2(77-77') , 2(77-77') 



sin[2/c(?7 — r/)] 
a 2 {rf)a 2 {rj') 
Id + Ind , 



k 3 r\ 2 r\' 2 

(vv') 2 

k^rfr]' 2 ^ k 2 m ' k 2 T] 2 f]' 2 



krjrj' 
2 



where and (j)t are the mode function and its complex 
conjugate respectively, given in Ea. ljAl|) . The only diver- 
gent term is 



)ad2 12 



1 



A 2 



87T 2 a 2 (77)' 

>iA ° (V) 



2^< 



2 -In 



a(r))n 



-In 



2A 



a(rf)fj, 



2A 



(B5a) 
(B5b) 



where 7 is the Euler's constant and we have replaced 
the bare parameters with their renormalized counterparts 



dk sin[2fc(?7 — 77')] 
fc>A 8vr 2 a 2 (77)a 2 (77') 



(B8) 



With the use of the definition of Id and Ind, we can 
write the order A 2 contribution as 



- 288 AV(t7) / drj'a^lMv) / d A y ^(i/)Jm(^ 
= 12A a 2 (77)0<(77) ([{<f>i)] D + m)Y ND ) , (B9) 
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with 

[^l)Y D = -24A J' 'do' a\o')4> 2 <{o')lD 

= J / <^ 27^ / dk sln 2fc W - »7 

J m a iV) Jk>A 

_ _ 3A 2 /" dfc 3A a 2 fa)0;lfa) 

- 2/< l,, L fc + 2** o»(»?) 

x / — cos[2fc(?7 - r?i)] 

3A /"' . (a^lfapy /■ dk 

+ ^ L " « a fr) Jk>A T cos[2fc(?? - " )] ' 



where the last equality follows after performing an in- 
tegration by parts. In this expression we can see that 
the first term after the last equality is (</>>) ^ aim hence 
all infinities arc cancelled. Finally, performing the Ind 
integral and reordering the terms we obtain the explicit 
form of the renormalized equation given in Ea. (|55|l . 
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